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Abstract 

A 3D singular integral equation is derived for electromagnetic wave 
scattering by bodies of arbitrary shape. Its numerical solution by a 
projection method is outlined. 
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1 Introduction 

Consider the following scattering problem. An incident electromagnetic field 
{Eq,Hq) is scattered by a bounded region D, filled with a material with 
parameters (e, a, /hq). The exterior region D' is a homogeneous region with 
parameters (eo,(T = 0,/xo). Consider for simplicity the case when e = const 
in D, a = const > in D. Let e' = e + ^. The governing equations in 
are 

V X E = iuJUoH, V X if = -iuoe'E. (1) 
At the boundary S of D one has 

[N,E+] = [N,E-], (2) 

and 

N -e'E^ = N -eoE-, (3) 

where is the unit normal to S, pointing into D', E^{E~) is the limiting 
value of on from inside (outside) S, [N^ E\ is the cross product, and 
E ■ N is the dot product of two vectors. 
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Let 



7 2 2 7^2 2 / v2 \ in Z)', , , 

fe = w eo/io, = a; e pto, /C = I ^2 ^ (4) 
Equations ([T]) imply 

V X V X - /c^^; = 0, zr = ^ ^ ^ m Ml (5) 

Therefore, it is sufficient to find E satisfying the first equation ([5]), boundary 
conditions ([2]), ([3]), and the radiation condition 

E = Eq + V] Vr-ikV = o , r := \x\ oo. (6) 

Equation ([5]) for E can be written as 

LE — VxVxE-k^E = pE; p ■= p(x) = JC^ -k^ = \ 

^ ' [ K"^ — k'^, m D. 

(7) 

The incident field Eq solves equation ([7]) with p = 0. 

( 1 i = j 

Let Six) denote the delta-function and Jj,- := < „' . , .' 

I 0, ^ / J. 

Let G = Gij{x) solve the problem: 

LG = 6{x)6ij , Gr-ikG = o , r ^ oo. (8) 

Then the solution to ([7])- ([6]) solves the integral equation 

E = Eo+ [ G{x-y)p{y)E{y)dy. (9) 

The kernel G{x) = G{\x\) is symmetric, Gij = Gji, see formula (fTSj) . Let us 
prove 

Lemma 1.1. There is at most one solution to Q satisfying Q and ([3]). 

Proof. If there are two solutions then their difference E solves the homoge- 
neous equation ([9]) and satisfies ([2]) and ([3]). Thus, E solves dl]), ([6]), ([2]) and 
([3]). Therefore, E and H = solve equations ([1]) and satisfy condition 
dl]), ([3]) and dSD. It is known (see e.g., [2j) that this implies E = H = 0. 
Lemma ll.ll is proved. □ 
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Lemma 1.2. If E solves equation then it satisfies ([2]), ([3]) and ([7]). 

Therefore, Q has at most one solution. 

Proof. Applying operator L to ([9|) one obtains equation ([7|). The integral 
in ^ is the term ^ in dH]). It satisfies the radiation condition because G 
does. Equation ^ is equivalent to ([5]). Equation ([5]) together with the 
formula H = ^^"^ yield both equations ([T]). Conditions ^ and (I3D are 
consequences of equations ([T|). Therefore, every solution to ([9]) is in one-to 
one correspondence with the solution to equations ([l]). This correspondence 
is given by the formulas E = E, H = By Lemma ll.ll equation ([9]) 

has at most one solution satisfying ([2]) and ([3]). We have proved that every 
solution to Q satisfies ([2]) and ([3]). Therefore, ([9]) has at most one solution. 
Lemma 11.21 is proved. □ 

Lemma 1.3. Equation ^ has a unique solution. 

Proof. Uniqueness of the solution to Q is proved in Lemma ll.2i Existence 
of it follows from the existence of the solution to the scattering problem and 
the fact, established in the proof of Lemma [1.2l that a solution to Q solves 
equation ([5]) and satisfies the radiation condition ([6]) and conditions ([2]), ([3]). 
Lemma 11.31 is proved. □ 

From lemmas 1-3 the following result follows 

Theorem 1.4. Equation ([9|) has a unique solution E. This solution E 
generates the solution to the scattering problem by the formula E = E, 



H = 



VxE 



In Section 2 we construct the Green's function G. 



2 Construction of G 



Let us look for G of the form 



G{x) = [ e- 



■^Gi^dC, G{0 



(27r)3 



1 




G{x)dx. 



(10) 



Take the Fourier transform of ([8|) and get 



-[e, [i^Gw-k'G 



1 




(11) 



(2vr) 
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where [a, b] is the cross product of two vectors, and a • 6 is their dot product. 
This imphes 

-^^■G+ie-k')G = j^I. (12) 



From (1121) one finds 



r< — (TA\ 



Thus, 

^ij ^i^j 

(2^)3(^2 _ (2^)3A;2(^2 _ ^2) 

Taking the inverse Fourier transform and using the radiation condition ([6]), 
one gets 

G^jix) = gix)6ij + -^dijoix); 5=^-^, di := —. (15) 
From p5|) and ([9]) one gets: 

^,(x) = Eo,{x) + (i^2 _ /■ ^(^^ y)E,{y)dy 

Jd 



— d f dg{x, y) 



(16) 



where summation over the repeated indices is understood. Equation (fT6|) is 
a vector singular integral equation. The operator 

TE={K^-k^) [ g{x,y)E{y)dy 
is compact in L^{D). Let 

^•=^^-^' QE = V I^V,g{x,y)-E{y)dy. (17) 



Then equation (jl6p can be written as 

E = Eq + TE + -iQE. (18) 

This is equation ([9]). 

Numerically one can solve equation (jl6p (or (jlSp ) by a projection method. 
For example, let {4>j{x)} be a basis of L'^{D) and € Hq{D), where Hq{D) 
is the closure of C^{D) functions in the norm of the Sobolev space H^{D). 
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Multiply equation (fT6]) by cp^ (the bar stands for the complex conjugate), 
integrate over D and then the third term by parts, to get: 



Eim — Eoim + {K — k 



M 

2 r,2\ 



/ / dx(p^{x)g{x,y) EimXpm'iv) 
Jd Jd 

dx^^^ I f2 E,m'cPm'{y)dy, 1 < m < M, 1 < ^ < 3. 

^ * J m'=l 

(19) 

This is a linear algebraic system for finding the coefficients: 

4^^^ ■■= E^m ■■= I E,{x)^^{x)dx. (20) 



JD 

The number M determines the accuracy of the appproximate solution E{x). 
One has 

hm \\E^'^^){x)-E{x)U^^j,)=Q. (21) 
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